Fractional and Integer Charges from Levinson's Theorem 



O 
O 



E. Farhi'', N. Graham'''^, R. L. Jaffe*^, and H. Weigel**^ 

^Center for Theoretical Physics, Laboratory for Nuclear Science and Department of Physics 
Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

and 

^Dragon Systems, Inc. Newton, MA 024-60 



MIT-CTP#3004 |hep-th/000718£ 



O 

o 
o 

^ Abstract 

Q 

We compute fractional and integer fermion quantum numbers of static 
CN ■ background field configurations using phase shifts and Levinson's theorem. 

. By extending fermionic scattering theory to arbitrary dimensions, we imple- 

OO \ ment dimensional regularization in a 1 + 1 dimensional gauge theory. We 



demonstrate that this regularization procedure automatically eliminates the 



' anomaly in the vector current that a naive regulator would produce. We also 



apply these techniques to bag models in one and three dimensions. 



■ . I. INTRODUCTION 

D , 

■ Many field theory solitons have especially interesting properties when they are coupled 

to fermions, because they act as strong background fields that can drastically alter the Dirac 
^ . spectrum. Solitons that break C and CP invariance can introduce asymmetries in the Dirac 
spectrum, causing the soliton to carry fermion number. By adiabatically turning on the 
soliton from a trivial background, one can observe this fermion number as a level crossing 
in the Dirac spectrum. In addition, solitons with nontrivial topological boundary conditions 
can carry fractional fermion number . 

Blankenbecler and Boyanovsky 0] showed how to use a phase shift representation of the 
density of states to calculate these fermion quantum numbers. There is no need to consider 
any interpolation of the background field. Together with the bound states, the behavior of 
the phase shift at threshold gives the integer charge. Fractional charges naturally appear in 
the behavior of the phase shift at large k. 

In recent work on variational computations of soliton energies |^^, it has been essential 
to know these fermion numbers, since one wants to compare configurations with the same 
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fermion charge. This work has led us to extend the method of Blankenbecler and Boyanovsky 
in two directions. First, we have shown how to extend dimensional regularization to a phase 
shift formalism. We are then able to reconsider the anomaly calculation of Ref. starting 
with an explicitly gauge-invariant regulator, rather than enforcing gauge invariance by hand 
on the result. We can thus explicitly see the difference between the naive regulator and 
the gauge-invariant regulator. Second, we demonstrate an extension of this work to chiral 
bag models in three spatial dimensions, giving a practical example of a fractional fermion 
number. 

Section II describes the Green's function formalism that relates the phase shifts to the 
Fock space expansion for the charge. This result can then be strengthened by the use of 
Levinson's theorem. In Section III we demonstrate the effects of the choice of regulator in 
a one- dimensional gauge theory. We show that by choosing a gauge-invariant regulator we 
automatically ensure that there are no anomalies in the vector current. In Section IV we 
perform explicit calculations in chiral bag models. In Appendix A, we develop the formalism 
for computing fermion phase shifts in arbitrary dimensions, and demonstrate that it leads 
to results that agree with the standard result from Feynman perturbation theory for the 
contribution of the tadpole graph to the effective energy. In Appendix B, we derive results 
used in Section IV for the three-dimensional bag. 



Our primary tool will be the scattering phase shifts. The phase shifts are useful because 
they unambiguously and quantitatively track the changes of the spectrum of small oscilla- 
tions around the background fields, even though that spectrum is continuous and infinite, 
and integrals over it of physical quantities may be divergent. 

To derive our results, we consider a single Dirac fermion in n space dimensions. We 
will assume that the soliton background has spherical symmetry, so that the spectrum 
decomposes into a sum over eigenchannels a. In one dimension there are just two channels, 
for even and odd parity, while in three dimensions the sum will run over parity and over 
either total spin or grand spin, also including the degeneracy within each channel. We begin 
by writing the fermion charge density as 



where x is an n-dimensional vector. We have been careful to order the anticommuting 
fermion fields to maintain charge conjugation invariance: is even under C and 75 ^E'] 

and [^,7''\E'] are odd under C. 

We make the usual Fock decomposition in terms of the eigenstates "0^ of the single- 
particle Dirac equation. 



II. SPECTRAL ANALYSIS IN SOLITON BACKGROUNDS 



A. Density of states 
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■) 



(2) 



—00 



2 



where the integral over lo includes both continuum and bound states. We have normalized 
the wavefunctions by 

y C(^)Va''(3;) dx = 6{cu - uj')6{a - a') (3) 

where the delta function of to on the right hand side is interpreted as a Dirac delta function 
for continuum states and a Kronecker delta for bound states. The creation and annihilation 
operators satisfy 

{K, b-^''] = {da, = ^(^ - - «') (4) 

with all other anticommutators vanishing. Thus the charge density becomes 

fi^) = -^E • (5) 

We next consider the Green's function for the fermion field 
G{x, y, t) = ^(n|T(v|/(x, t)^{y, O^m 

= ^ E (/ rf^e-*C(x)C(2/)^e(-t) - da;e— (6) 

and its Fourier transform 

whose trace gives the density of states according to 

p(oo) = sgn(a;) = Im Tr - / G{x, x, co) (8) 

giving as a result 

'^H = ^E / ^"^ic(^)r- (9) 

a 

for > m. (For |a;| < m, the density of states consists of a sum of delta functions for 
each bound state.) A more useful form of this equation is obtained by subtracting the free 
density of states p°(a;), giving 

- p'iu;) = ^ E / ^"^(IC(^)r - 1) ■ (10) 

a 

We can also express the density of states in each channel in terms of the S-matrix: 

/\ n/\ / \ ^ d r \ I \^ d^cSj^) /-.-.\ 

PocK^) - PcM) = ^S^(^)^^^°S'^"('^) = sgn(cc;)- . (11) 
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The total density of states is then 

p{lo) = Y,Pc.{^)- (12) 

Putting these results into eq. we thus obtain an integral over the continuum and a 
sum over bound states 

duj dSai^jj) 1 duj d5a{u>) sr-^ 1 | 

2tx duj ^ ^ 2 ~ 27r du ^ 2 I 

= — (5„(m) - 6a{oo) - Tvn^ + rrn^ - ^^(-m) + 6a{-oo)) (13) 
Zn ^ — ^ 

a 

where and n< give the number of bound states with positive and negative energy respec- 
tively in each channel. This is our main result. We will see that it generalizes unchanged to 
cases with fractional charges, which appear in the phase shift at u; = ±oo. 




B. Levinson's Theorem 

We have seen that the phase shifts keep track of the rearrangement of the fermion 
spectrum in a quantitative way, by telling us the density of states. Levinson's theorem 
shows that the phase shifts also track how many states have left the continuum at each 
threshold. This property will give us a physical motivation for eq. (p!3[), and will also enable 
us to simplify it. 

Levinson's theorem in three dimensions and in the negative parity channel in one di- 
mension gives the number of states that have left the continuum by passing through the 
threshold at m as 

6{m) - 6{oo) = Nn. (14) 

These states typically appear as bound states (which give delta functions in the density of 
states), though it is possible that in cases where the spectrum is not charge conjugation 
invariant, they can reenter the continuum of states with opposite energy. In the positive 
parity channel in one dimension, Levinson's theorem is modified to 

5(m) - (5(oo) = (TV - i)7r. (15) 

In one dimension and in the i = channel in three dimensions, there is the possibility of a 
state whose wavefunction goes to a constant at r = oo, rather than decaying exponentially 
like a bound state or oscillating like a continuum state. Such "half-bound" states are on 
the verge of becoming bound and count as | in Levinson's theorem. In the free case, for 
example, S{uj) = for all u, but eq. (p!5D still holds because there is a half-bound state with 
ijj = constant. For fuller discussion of these results, see Ref. 0. 

Computing the fermion number of a field configuration now becomes a matter of simple 
counting. We consider each channel separately. If a bound state leaves the positive contin- 
uum but appears as a positive energy bound state, it has not changed the fermion number 
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of the configuration. However, if it crosses 00 = and becomes a negative energy bound 
state, it is now filled in the vacuum and gives a fermion number of one. Thus the fermion 
number of a field configuration is given by 

Q = 77- i^cim) - 6a{oo) - 7rn> + 7rn< - (5„(-m) + 6a{-oo)) (16) 
Ztt ^ — ' 

a 

which is exactly (plsl) . 

This interpretation suggests a further simplification: since Levinson's theorem tracks all 
states that enter and leave the two continua, even in the presence of CP-violation, we have 
the restriction that 

6a{m) - 5a{oo) + 5a{-m) - 5a{-oo) - nn^ - Tvn^ = , (17) 

so that 

Q = \^ - {5aim) - Saioo) - 7rn>) 

^-^ 71 

a 

= V-(7rn<-5,(-m) + <5„(-oo)). (18) 



TT 

a 



In the positive parity channel in one space dimension, we must subtract 1 from the left- 
hand side of eq. ([TtD and | from the subsequent expressions for Q because of the modification 
to Levinson's theorem in eq. (0). 



III. ELECTROSTATICS AND THE NEED FOR REGULARIZATION 

The conserved charges we consider are not renormalized. That is, they do not receive 
any contributions from the counterterms of the theory. Nonetheless, it is essential to include 
the effects of the regularization procedure used to define the theory. The example of QED in 
H-l dimensions provides a clear illustration of this subtlety. Although the theory is finite, 
the regularization process is nontrivial. 

The Lagrangian is 

C = -^F.^F'^" + \ {^{^9, - A,) - m) ^] (19) 

where we again have used the commutator to ensure that the free theory is C and CP 
invariant. In ordinary perturbation theory, the vacuum polarization diagram computed in 
d space-time dimensions is 

n r ^ 9- f\,f 2C(1 - Oig,.p' - p.p.) + g,. jm' - p'gi - Q + P(| - 1)) 
n,.(,) = 2^e N,j^ d^j ^ (F+^^^d-O-^T 

(20) 

where 2Na is the dimension of the Dirac algebra. If we had not regulated the theory by 
analytically continuing the space-time dimension, we would not have found the last term. 
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which vanishes if we set d = 2 from the outset. Keeping d ^ 2 shows that this term exactly 
cancels the two terms that precede it, leaving the transverse form of the vacuum polarization 
that is required by gauge invariance. Thus we must include in our definition of the field 
theory the additional information that the theory is regulated in order to preserve gauge 
invariance at the quantum level, and dimensional regularization provides a convenient way 
to implement this requirement. 

The vacuum polarization diagram reflects the effect of the anomaly. The anomaly is 
obtained from the leading correction to the vector current, which is related to the polarization 
tensor by 



j,{x) = / An^.(x - y)A''{y) (21) 



where Ii^y{x) denotes the Fourier transform of eq. (^). Thus a completely transverse 
polarization tensor corresponds to a conserved vector current. Setting d = 2 from the outset 
gives an anomalous vector current 

d,j^ = a^^7^75^ =^0 . (22) 

Including the contribution proportional to | — 1 in eq. (|20| ) transfers this anomaly to the 
axial current 

d^f = d^^Y^ = 

d,j^ = d,^Yl5^ = ^'"'F.u . (23) 

If we choose a configuration with Ai = and adiabatically turn on a configuration Aq{x) 
between t = — oo and t = 0, then integrating eq. (|2^) we obtain 

Q = J f{x)dx = ^ J Mx)dx (24) 

for the naive regulator while Q = in dimensional regularization. 

The phase shift approach shows exactly the same behavior. We consider the example of 
an electrostatic square well potential with depth tp and width 2L. First we ignore subtleties 
of regularization and compute directly in c? = 2. The phase shift in the negative parity 
channel S~{uj) is determined by 

— - — tanikL + 5" (iu)) = -tanqL (25) 

k q 

with k = y/u"^ — m? and q = ^J{lo + eipY — m"^- Similarly, the phase shift in the positive 
parity channel S~^{uj) is determined by 

^ tan(A;L + 5+ {to)) = tan qL . (26) 



m + u m + u + eip 

As u ^ ±oo, the total phase shift approaches ±2eLpL, giving a fractional contribution to 
the total fermion charge in agreement with eq. (p4|). Although such fractional charges are 
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possible (and we will see examples of them in the next section), in this case the result 
indicates that the method of calculation has not preserved gauge invariance. 

To preserve gauge invariance, we regularize the computation by computing the phase 
shifts as analytic functions of the space dimension n. Since we are only concerned with 
the contribution from u —>■ ±00, we can consider just the leading Born approximation. In 
Appendix A we have extended the method of dimensional regularization of the phase shifts 
in Ref. to fermions. In arbitrary dimensions, the total phase shift is a sum over channels 
labeled by total spin j = |, | . . . and by parity. Summing over parity, the leading Born 
approximation to total phase shift in each j channel is given by 

roc 

^2 = uen / Ao(r) (j.+,_3(fcr)2 + J„_,^_i (/cr)^] rdr (27) 



which has degeneracy d{j) given by eq. ( |A19| ). Summing over j using eq. ( |A16| ) and eq. ( |A23[ ) 



yields the leading Born approximation to the total phase shift in n space dimensions 



roc 

ue-KNdY.W) / Ao(r)J|+,_i(A;r)Vt/r 



- 2"-2nr(f )2 

which reduces to N^eifL if we send n 1 and take the limit as a; — ±00. But the order 
of these limits is essential: if we first regulate the theory by holding the dimension fixed at 
ri < 1, we then see that the contribution as a; — ±00 vanishes. Only after we have taken the 
uj ±00 limits do we send n ^ 1. This procedure, dictated by dimensional regularization, 
preserves gauge invariance and gives no fractional charge. 

We note that other regularization methods commonly used in phase shift calculations, 
such as zeta-function regularization, would not preserve gauge invariance and would thus 
lead to the same spurious fractional result. 



IV. FRACTIONAL CHARGES 

A. Chiral bag model in one dimension 

Chiral bag models provide simple illustrations of fractional fermion numbers. We begin 
with a Dirac fermion in one dimension on the half-line x > 0, subject to the boundary 
condition 

^gnss^ = (29) 
at X = 0, with — I < 9 < |. We consider the Dirac Hamiltonian 

7°(-i7ia^ + m)m = io^ (30) 
in the basis 7° = and 7^ = ia2. The free solutions are 
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0±(x) = ( ) e^*'^" (31) 



uj — m 



where k = y/u"^ — m?. The full solutions are then given in terms of the phase shifts by 

^± (x) = 0- (x) + e^'^^ (2;) • (32) 
We can then solve for the phase shifts using the boundary condition, eq. (^91). We obtain 



k 

cot5"'"(a') = tan/? 

LO — m 
k 

tan5^(tt;) = tan/? (33) 

u — m 

where /? = f — To find the bound states, we look for solutions of the form 

m={ (34) 



with K = \Jrn? — u^. Imposing the boundary condition gives 

K = (m + cj) tan/3 (35) 

so there is always exactly one bound states for —\<Q<\- Plugging these results into 
eq. ([I3D , we find that the fermion number is ^. 

B. Chiral bag model in 3 dimension 

This simple model generalizes naturally to 3 dimensions. We consider an isodoublet of 
Dirac fermions subject to the boundary condition 

^e^^"-"^^^ = (36) 

imposed on a sphere of radius i?, where r are the isospin Pauli matrices. This condition 
is not invariant under space and isospin rotations individually, but it is invariant under 
combined space and isospin rotations, and under parity. Thus we can decompose the scat- 
tering problem into eigenchannels labelled by grand spin G = 0, 1, 2 . . . and parity. This 
calculation is outlined in Appendix B. 

In the G = channel, the phase shifts in the two parity channels are 

btiuj) =ATg(ihi(kR) cose + tho(kR)(lTfiinO)] (37) 

V ±00 — m I 



where hn{x) are spherical Hankel functions of the first kind and k = \/ uP' — . Extracting 
the contribution from the G = channel is straightforward. Using eq. ([13|) , we find a 
contribution of — - to the fermion number. 

For G > there are two states for each choice of G and parity. The total phase shift in 
each channel is given by 
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5^{uj) = Arg| sin^ .^^^^^.^' + ^ ^ hcikR) {hc+iikR) - hc-iikR)) 
' [uj^m)R ±.uj — m 



- cos' d 



hG+i{kR)hG-i{kR) 



k 



±0; — m 



hcikRf 



(38) 



The contribution from G > needs to be treated with care, 
contribution comes from G + ^ 



For large k, the dominant 
2 - kR. We must follow a consistent regularization procedure 
in order to obtain the correct order of limits. In this model a simple cutoff suffices. We 
first compute the total phase shift at = A, summed over all partial waves. This sum then 
has a smooth limit as A ^ oo. If we had taken the limit in the other order by considering 
the A; — >■ oo limit in each partial wave separately, we would incorrectly conclude that the 
contribution from G > was identically zero. 

For numerical computations, it is convenient to first consider ^ and then integrate to 
obtain Q as a function of 9. Figure |^ shows the ^ computation for one value of 6. Fixing 
a large value of the cutoff A, we then sum over both parities from G = up to Gmax ^ Ai?, 
so that the contribution from the higher values of G is negligible. We then compute the ^ 
for each 9 from eq. (O) using A; = A in place of A; = oo. 




FIG. 1. Numerical computation of ^ J2g='o('^G + 1) {6^{io) + S^{uj) - S^{-u) - 6^{-uj)), 
evaluated at 9 = ^ and m = 0, summed up to Gmax = 32, and plotted a function of z = ojR for 
a; > 0. The larger values of G that we have omitted are negligible for z <C Gmax- Tor z « Gmax/2 
we see good agreement with the result of Ref. which is shown hy a horizontal line. For z ^ Gmax 
we see only the contribution from G = 0, which would lead to an incorrect result for the charge. 

Carrying out this computation numerically for each 9, we find agreement with the es- 
tablished result 

§^-i(l-co=2.,. (39) 

SO that 

Q = --(^0 -sm9cos9). (40) 
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V. CONCLUSIONS 



By explicitly introducing appropriate regulators, we have been able to unambiguously 
compute fermion quantum numbers in a variety of models. We have seen that by introducing 
the regulator as part of the definition of the theory, we have no need (and no freedom) to ad- 
just the results of our computations after the fact. This determinism leads to a computation 
that is both theoretically elegant and computationally practical. 
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APPENDIX A: DIRAC EQUATION IN ARBITRARY DIMENSIONS 

In this Appendix we discuss the Dirac equation for a radially symmetric potential 
V — V{r) for an arbitrary number of spatial dimensions n. We calculate the leading Born 
approximation to the phase shift in n dimensions and find its contribution to the Casimir 

energy. We verify that this result agrees with what we obtain from the tadpole graph in 
standard Fcynman perturbation theory using dimensional regularization. 
Our starting point is the free Dirac equation in n spatial dimensions, 

(^-ia ■ d + (3rnj * = = u;* . (Al) 

The spinor ^ has 2Nii components and accordingly the Dirac matrices aj,j = 1, . . . ,n and 
p have 2Nd x 2Nd elements, where = 2^''-^^^ for n odd and = 2"/^ for n even. We 
will generalize the case of n odd, though our results will not depend on this choice. We 
choose the basis 

and J^l,...,n. (A2) 

The Clifford algebra is obtained by demanding the anti-commutator {Aj,Aj} = 25ij. In 
analogy to the spin generators in n = 3 we define the commutator 

[Ai,A,-] = 2iEi,- (A3) 

which obeys the SO{n) commutation relation 

[Ejj, Efe;] = i {Sik'Sji + SjiTtik — Sii'Ejk — 5jk^ii) . (A4) 

We define the orbital angular momentum operator 

Lij = -i {xidj - Xjdi) (A5) 
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which also satisfies the SO{n) algebra 



[Lij, Li:i] — i (SikLji + 5jiLik - SiiLjk — SjkLi 



(A6) 



We can then put these together to form the total spin operator 



(A7) 



which commutes with the Hamiltonian: [H, Jij] = 0. Having obtained the algebra, we next 
need to find the Casimir eigenvalues of 



and 



9 5Z "^S • 



(A8) 



The eigenvalues of are those of SO{n), i{i + n — 2). To find S^, we consider its trace, 
which is just the number of independent matrices Sjj, ^n{n — l). Then to obtain the Casimir 
eigenvalue of J^, all we need to find is {L ■ T,) = j LijUij) . We use the second order 

equations obtained from eq. ([Al|) , which are generalized Bessel equations. We first remark 
that f ■ A has zero total spin ([f • A, Jjj] = 0) and eigenvalue i = 1 with respect to L^. 
Therefore the appropriate spinor with definite parity can be parameterized as 



g{r) (r ■ A) 1^,,,,- 



(A9) 



where l^,s,j denote generalized spinor spherical harmonics. The radial functions obey the 
coupled first order equations 



n-l + R 
2r 

1 - R' 



n 



2r 



g(r) = (m - uj)f{r) 
fir) = {m + u)g{r) 



(AlO) 



where R = (L ■ S) + n — 1 contains the desired eigenvalue. We can decouple these equations 
to obtain second order equations for /(r) and g{r). By demanding that /(r) and g{r) obey 
generalized Bessel equations with orbital angular momentum i and £' respectively, we find 



R = l±{n-2-2i) 



and 



R 



-1 ± (n - 2 + 2£') . 



(All) 



In view of the above mentioned properties of f - A we have i' = £±1. Hence the two relations 
in eq. ( [All|) are consistent if 



R = n + 2i-l for £' = £+! and R = 3 - n ~ 2£ for 
Putting these results together we find the Casimir eigenvalue 



- 1 . 



= i(^i + n-2) + 



n{n — 1) 



£, £' = £+! 

2-n-£, £' = £-! 



(A12) 



(A13) 
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Defining j = h{i + i') yields 



J^=l,-iU,+„-!U!l(^ (A14) 



2/ V 2^ 

for botli cases. Tlie above definition of j also ensures that (as for n = 3) there are two 
independent solutions for a given j: {}) ^ = j + \A' = 3 ~ \ and (ii) d- = j — \A' = i + \- 
These solutions have opposite parity. 

Finally, we have to find the degeneracy factor, d{j). We use two trace relations that are 
simple in the basis appropriate for and Lij. Written in the basis for Jij they connect 
different representations and provide information about the degeneracy factors. We have 

d{j){R{3)-n + l) = Q (A15) 



and 

dU) = N,D{i) . (A16) 

where we have defined c?(— |) = independent of n. The first condition is nothing but the 
tracelessness of {L ■ S) while the second gives the number of states for a given orbital angular 
momentum i, 



Eq. (|A15|) can be re-expressed as 



dii+-) = ^^^^d{i-^) (A18) 
which after substitution into eq. ( [A 161 ) yields the final result 

d(j) = NAj + -) , \ , . (A19) 

Equation (|A1^ ) represents a recursion relation between di^f) and di^j + 1) that can straight- 
forwardly be shown to be satisfied by the degeneracy factor of eq. (|A19|) . We note that 
as n — > 1, eq. (|A19|) gives zero in all channels except 3 = \i where it is one. As in the 



bosonic case, this limit gives the reduction to the positive and negative parity channels in 
one dimension. 

To show how the dimensional regularization of the phase shifts corresponds to ordinary 
dimensional regularization of Feynman perturbation theory, we compute the contribution 
to the Casimir energy from the leading Born approximation to the phase shifts using the 
method of |Q, and compare it to the energy of the tadpole graph, with both quantities 
evaluated in n dimensions. 
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The leading Born approximation to the phase shift is, summing over parity channels, 

= -2 / d^^^^y [J^+j-likrr + J.^j-iikrf) (A20) 

«/ 

for k — \/uP- — m?. The leading Born approximation to the Casimir energy is then given by 

1 poo (^(^(1) 

A£;(i)[0] = — / dkiVk-' + m2 - m) V d{j)-^ . (A21) 
TT 7o J dk 

Summing over all channels we find the total phase shift 



J 2 ' 2 • 



--irdrV{r)r ^ d{j) (ju^,_.{krf Jn^,_.{krf) 

Jo 1 3 

J 2'2'" 

/'OO / -| -| 

= --/ dry(r)r Ui + -)J.+e-i{krf + d{i - -)Jn^e-i{krf 

•^^ £=0,1,2... ^ 

/"OO 

= --/ drV{r)r ^ 7V,D(£) Jn+,_i(A;r)2 . (A22) 

£=0,1,2... 



Using the Bessel function identity 

(2g + 2£)r(2g + £) 2 r(2g + l)/z\29 



^ fxZTT) ^rx^TTF^sJ ^"^^^^ 

and setting g = | — 1, we sum over £ and obtain for the Casimir energy 

A£;(i)[^] = 2N, ^Yl.A n - 2) H {uj - m)k^-^ dk (A24) 
(47r)2r(f) Jo 



where 



r 27rt /"°° 

{V) = / = —r^ / y(r)r"-^dr . (A25) 

The /c integral can be calculated in the vicinity of n = | and then analytically continued, 
yielding 

f(.-™).-3,,._5,(i^),(^). ,,,, 



Hence we find 



A£;« [0 ] = -2Ar,^^r (^—^) m--' (A27) 
(47r)~ \ ^ J 

which is exactly what we obtain using standard dimensional regularization for the tadpole 
diagram in n + 1 space-time dimensions. 
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APPENDIX B: CHIRAL BAG MODEL S-MATRIX 



In this Appendix we outline the derivation of eq. (P7|) and eq. (PSp. We begin with 



spinors that are eigenstates of parity and total grand spin where grand spin G is the 
sum of total spin j = l + and isospin ^r. For a given grand spin G with ^-component M, 
we will find the scattering wavefunctions in terms of the spherical harmonic functions yj^i 
with J = G ± ^ and £ = j ± ^. These are two-component spinors in both spin and isospin 



space. While grand spin is conserved, the boundary condition in eq. (|36D mixes states with 
different ordinary spin j. For the channels with parity {—)^ we have two spinors that solve 
the Dirac equation away from the boundary: 

For zero grand spin the second spinor is absent. We now introduce linear combinations that 
define the S-matrix: 



w-hG-i{kr)ya_iG-J V-^"^G-i(^03^G-iG-i 



iw~^hr(kr)yr,ir \ , . 

\w h*a^^{kr)ya+i^a+i J 

where hi{kr) refers to the spherical Hankel functions suitable to parameterize an in- 
coming spherical wave. We have introduced the kinematic factors = a/^H-^ and 
w~ = sgn(a;)iyi — ^, where u = ±yk^+m? and m are the energy and mass of the Dirac 
particle respectively. In the case G = 0, the components with j = G — \ are absent and the 
S-matrix has only a single component 5*^'° = exp(2i5(]"). 

Imposing the boundary condition, eq. (p6D, on these wavefunctions gives 



cos^ ir ■ fsinO — if ■ a 

if ■ fsixiO + if ■ a cos6 



n 

r=R 



for i = l,2. (B4) 



For each n = 1,2 the projection onto grand spin spherical harmonics yields two equations, 
which allows us to extract all four components of the S-matrix. It is convenient to express 
the result in the form of a matrix equation: 

where the star denotes complex conjugation. The components of the matrix M are given by 
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X = h*(kR)cose + — ^ (1 + — ^ — sin^ I h* (kR) 
CO + m \ 2G + 1 / 

X = -ha{kR)cose + — ^ (1 + J' sine] hG_^{kR) 
2v/G(GTi) k 

The total phase shift in this channel is then given by 

6% = Arg (det^-^'^) = Arg (B7) 

leading to eq (|38|). For the channels with parity {—)'~^^^ the computation proceeds analo- 
gously from the definition of the corresponding S-matrix: 

^ / iu;+/iG+i(/i;r)3^G+i,G+i\ _^^-,G f iw+h*G^^{kr)yc+^ c+i 



-w hG{kr)y^_^iQ J \ -w hUkr)yG+^^G 



+5.7"f^"''i;i^\^"--"-M (B8) 



iuj+hG-iikr)ya_i^G-i \ _^ ^-,g ( ^«'^^G-i(^03^G-i,G-i 

w'hGyG-^n J V w-h}.y^_.(. 
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